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Abstract 

q Quantum field theory violates all the classical energy conditions of general rela- 

tivity. Nonetheless, it turns out that quantum field theories satisfy remnants of the 
classical energy conditions, known as Quantum Energy Inequalities (QEIs), that have 
been developed by various authors since the original pioneering work of Ford in 1978. 
These notes provide an introduction to QEIs and also to some of the techniques of 
quantum field theory in curved spacetime (particularly, the use of microlocal analy- 

<3\ sis together with the algebraic formulation of QFT) that enable rigorous and general 

QEIs to be derived. Specific examples are computed for the free scalar field and their 

l^-j consequences are discussed. QEIs are also derived for the class of unitary, positive 

energy conformal field theories in two spacetime dimensions. In that setting it is 
also possible to determine the probability distribution for individual measurements 
of certain smearings of the stress-energy tensor in the vacuum state. 
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Introduction and scope 



The weak energy condition (WEC) of classical general relativity holds if all observers 
measure the local energy density of matter to be nonnegative. Expressed in terms of the 
stress-energy tensor T a b, the WEC amounts to the inequality T ab u a u h > for all timelike 
vectors u a . This condition, and others like it, were introduced as conditions to model 
generic matter distributions in the discussion of results such as singularity theorems. The 
energy conditions are obeyed by many (though not all) classical theories of interest and 
have a central place in mathematical general relativity. 

However, it has long been known that quantum field theory, our most successful fun- 
damental matter model, is incompatible with these energy conditions [14]. This raises 
many questions, for example: Should one doubt the validity of the singularity theorems 
for realistic matter? Can quantum fields be used to support 'exotic' spacetime geome- 
tries (time machines, wormholes, warp drives...) which require energy-condition-violating 
stress-energy tensors if they are to be solutions to the Einstein equations? 

These lecture notes provide an introduction to the subject of Quantum Energy In- 
equalities (QEIs) [also often called Quantum Inequalities (QIs) in the literature]. These 
are conditions derived from within quantum field theory, that constrain the extent to which 
the classical energy conditions are violated. They strongly suggest that the answers to the 
questions just asked are negative (or raise apparently infeasible engineering problems). 

The lectures were given at the Albert Einstein Institute, Golm (March 2012) as five 
lectures aimed at graduate students enrolled in the IMPRS programme. The secondary 
aim of the lectures was to provide an introduction to the algebraic formulation of quantum 
field theory in curved spacetimes and the microlocal analytic techniques that, following 
their introduction by Radzikowski |74j . have been important in many recent developments, 
including the perturbative construction of interacting field theories in curved spacetime 
[58j EH [56]. From this perspective, the QEIs provide a useful application of the theory 
that illustrates some of its key features. With this in mind, the discussion of QEIs is 
biased towards the derivation and analysis of various bounds, and not to applications 
such as constraints on exotic spacetime geometries. The reviews [751 EH H2] provide a 
counterpoint in that regard. 

I would like to thank the organisers of the IMPRS lectures, particularly Stefan Fre- 
denhagen, for the invitation and for financial support, the students for their interest and 
questions, Atsushi Higuchi and Henning Bostelmann for useful discussions and Tom Roman 
for comments on the text. 

Summary of main conventions 

• h = c = G = 1, and the metric has signature H 

• (V a V 6 - V 6 V> d = RahM, and hence (V a V b - V h V a )v d = R abcd v c ; 

• the Ricci tensor is R ab = R d a db 

• Fourier transforms will be defined nonstandardly by 




the hat will sometimes be displaced e.g., f A {k), for typographical reasons. 
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1 Quantum (energy) inequalities 



1.1 The classical energy conditions 

The Einstein equations 

Gab = Rab — -jRQ 'ab = ~ 87rT a b 

are a curious mixture. Einstein himself likened the theory to a building 'one wing of which 
is built of fine marble... but the other wing of which is built of low-grade wood'. Elsewhere, 
he wrote that 'the right side is a formal condensation of all things whose comprehension 
in the sense of field-theory is still problematic' [j] 

Two observations illustrate some of the 'low grade' features of the right-hand side. First, 
non-gravitational physics only ever considers differences between stress-energy tensors (or 
derivatives of stress-energy tensors). Arguably, the quantity on the right-hand side of the 
Einstein equations [the absolute stress-energy tensor] appears nowhere else in physics^] 

Second, without further qualification, the Einstein equations have no predictive power 
whatsoever: every smooth Lorentzian spacetime solves the Einstein equations for a suitable 
choice of T ab - a fact that is often exploited in discussions of time-machines, wormholes, 
warp drives etc. A key issue is to determine what sorts of stress-energy tensors are physi- 
cally reasonable. As the real world contains a complicated combination of many different 
sources of stress-energy, one would ideally like some general principles that should apply 
fairly widely. 

The energy conditions are attempts at such principles. They are motivated partly by 
physics, and partly by mathematical expediency. The main conditions are: 

• The Weak Energy Condition (WEC) 

T ab u a u b > for all timelike u a 

Interpretation: All observers see nonnegative energy density. 

• The Null Energy Condition (NEC) 

T ab u a u b > for all null u a 

• The Dominant Energy Condition (DEC) 

T ab u a v b > for all future-pointing timelike u a and v b 

Interpretation: All observers see a causal flux of energy-momentum. 

• The Strong Energy Condition (SEC) 



T ab u a u b - l -g ab T ah > 



for all timelike unit u a . 



Quoted in [TU] Ch. 5, p. 123. 

2 Thc standard approach-which we employ later-is to obtain the stress-energy tensor as a functional 
derivative of the action with respect to the metric. But many actions describe the same physics, of course. 
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If matter satisfies these conditions and the Einstein equations hold, we obtain corre- 
sponding conditions on the geometry. For example, the NEC implies that R ab u a u b < for 
all null u a , while the SEC implies the same for all timelike u a . 

The principal interest in these conditions comes from the fact that they enforce various 
focussing behaviours for congruences of geodesies (see, e.g., [54, 86J). Consider a congruence 
of timelike geodesies with future-pointing unit tangent field u a . The expansion, shear and 
vorticity are defined uniquely by 

V b n a = ^9(u a u b - g ab ) + o ab + u ab 

and the requirements that a is traceless symmetric, while u is antisymmetric. In particular, 
= V a u a , and a ab is purely spatial and has positive square a ab a ab . The Raychaudhuri 
equation)^] gives 

^ = R ab u a u b - \e 2 - a ab a ab + co ab co ab 
dr 6 

so it is clear that an irrotational congruence obeys 

d9 1 



if the SEC holds, or equivalently that 



Thus 



and hence, 



dr 3 



d 1 1 

dr6 ~ 3 



1 1 T 

> 



6(t) ~ 6(0) 3 
(t) < ' 



fl(O)- 1 +T/3' 

Accordingly, if 0(0) < 0, we see that 0(r) -> — oo as r -> r t < 3|0(O)| _1 . Initially 
contracting irrotational congruences therefore form focal points in finite proper time if 
the ambient matter distribution obeys SEC. Arguments of this type lie at the core of the 
singularity theorems [51] and many other key results in mathematical relativity. 



3 Derivation: Use the geodesic property of u and the Riemann tensor definition to obtain 

u c V c V 6 M a = w c V 6 V c u Q - R cba d u c u d = V 6 (w c V c u Q ) - (Vf,u c )(V c w a ) - R cba d u c u d 
= -(V h u c )(V c u a ) - R cba d u c u d 

and trace to give the required result, noting that the terms in the expansion of Vbii a are mutually orthog- 
onal. 
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Examples 

1. A perfect fluid has stress-energy tensor 

Tab = (A* + P)^6 - PfiW, 

for flow 4-velocity t> a , energy-density fi, pressure p. Noting that the contraction u a v a 
between future-pointing timelike unit vectors u a and v a may take any value in [1, oo), 
one may easily show that 

• WEC holds iff fj, > and \i + p > 

• SEC holds iff fj, + p > and fi + 3p > 

• NEC holds iff fi + p > 

• DEC holds iff > \p\. 

Exercise: Prove these statements. Deduce that, while it is clear that DEC 
WEC =^ NEC, and also that SEC =^ NEC, there are no implications between 
SEC and DEC/WEC. 

2. The minimally coupled scalar field has stress-energy tensor 

T ab = (V a 0)(V 6 0) - ^a^ cd (V c 0)(V d 0) + ^mW 2 
Given any timelike unit vector u a , choose a tetrad with eo = u. Then 

3 

T afeM V = - J2( e A ■ V0) 2 + -m 2 2 > 

so this theory obeys WEC due to the 'sum of squares' form. Likewise, DEC also 
holds for this reason, and the same is true for the electromagnetic field. 

However, we also see that 

T ab u a u b - l -T\ = {u a V a ^) 2 - \m 2 tf 

so the SEC can fail even for this model if m > 0. 

3. The nonminimally coupled field, with coupling £, has stress-energy tensor 

TS = T { a t° ] + £ (<**□, - V a V b - Gab) 4>\ 

As the additional terms are not of the sum of squares form, even NEC can be violated: 
at points where V0 = 0, for example, this happens when the second derivative terms 
outweigh contributions proportional to 2 . 
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1.2 Violation of Energy Conditions in QFT 

Quantization and positivity do not mix well. For instance, the prototypical example of a 
positive classical quantity — the square of a field — is replaced by a Wick square in QFT. 
Although the formal square is positive, it is infinite; although the Wick square is finite, it 
is indefinite. 

To be specific, consider the standard quantized real scalar field in four-dimensional 
Minkowski space with mass m > 0, 



r rl 3 k 

$(*) = J (e-^a(fc) + e*-a(fc)*) 



(27r)V2w 

where the 4-vector k is k^ = (uj,k), with oj = ^/\\k\\ 2 + m 2 , and the annihilation and 
creation operators obey the commutation relations 

[a{k), a{k')] = 0, [a(k), a(k')} = (2vr) 3 5 (3) (fc - k')t. 

The vacuum vector, annihilated by all the a(k), is denoted Q. To form the Wick square 
:$ 2 :(a;), of course, we modify the formal expression for $(a;) 2 , replacing a(k)a(k')* by 
a(k')*a(k). We define the smeared Wick square by 



:$ 2 :(/) = J r/ ; ,-:<I. 2 :(,-)/i.r). 



where / G C^°(1R 4 ) is any test function on spacetime. Then it is a simple calculation 
(recalling our convention for Fourier transform) to show that 

J [Ik) 6 (2tc) 6 2Vuju' 
it is obvious that (Q | :$ 2 :(/)0) = 0, and a short calculation gives 

lid. 2 moil 2 = f d " k d ' k ' \f( k + k ')\ 2 

11 - {J) 11 7 (27r) 3 (2vr) 3 2uu> ' 

which is nonzero unless / is identically zero^J The observable :$ 2 :(/) therefore has van- 
ishing expectation value in the state Q, but does not annihilate Q. Standard variational 
arguments imply that :$ 2 :(/) must have some negative spectrum. Indeed, if we write 

ip a = cos a f2 + sin <x$ 2 : (/)f2 

(assuming / is chosen so ||:$ 2 :(/)0|| = 1) it is easy to calculate 

(V> I :$ 2 :(/)^a) = sin(2a) + sin 2 a(tt | :$ 2 :(/) 3 0) = 2a + 0(a 2 ) 



^Indeed, this is true on general grounds owing to the Reeh-Schlieder theorem. 
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giving negative expectation values for sufficiently small a < 0, even if / is nonnegative. 
By a scaling argument [18] it may be shown that the expectation value of :$ 2 : at a point 
is unbounded from below as the state varies among Hadamard states. 

The same conclusion may be reached on general grounds. An argument due to Epstein, 
Glaser and Jaffe [Hj proves that loss of positivity is unavoidable for Wightman fields with 
vanishing vacuum expectation values. The main thrust of their argument is the following. 
Suppose a local observable A has vanishing vacuum expectation value, i.e., (Q \ AQ) = 0. 
If A is positive, it has a square root, and we have 

p 1 / 2 ^!! 2 = (Q | AQ) = 

and therefore A l / 2 Q = 0. Hence AQ = and, as the Reeh-Schlieder theorem (52] tells us 
that no nonzero local observable can annihilate the vacuum, A must vanish. (The Reeh- 
Schlieder theorem only applies to local observables, which is why there is no contradiction 
between the positivity of the Hamiltonian and its vanishing v.e.v.) Alternatively, one can 
argue as follows: individual measurements of A in state Q constitute a random variable 
with vanishing expectation value; this implies either that zero is measured with probability 
1, in which case AQ = (impossible for nonzero local observables by the Reeh-Schlieder 
theorem) or that there is a nonzero probability for both positive and negative measurement 
values, so the spectrum of A extends into the negative half-line. 

There are many physical situations of interest in which negative energy densities arise 
in QFT calculations. One of the main examples is provided by the Casimir effect, in which 
plane parallel conducting plates in vacuo experience an attractive force. Actually, quite a 
bit can be done without much calculation [6] . In the case of infinite plane plates, separated 
through distance L along the z-axis in standard (£, x, y, z) coordinates, one may deduce 
on symmetry grounds and dimensional considerations that the stress-energy tensor of the 
electromagnetic field takes the form [6] 

T M , = ^diag (-1,1, 1,-3). 

where C(z) is dimensionless; here we have also used tracelessness of the stress-energy 
tensor. Conservation of the stress-energy tensor entails that C(z) is constant except at the 
plates, so C may take different values C and C\ inside and outside the plates (by reflection 
symmetry the values on the two outer components should be equal). As there is no other 
length scale in the problem, C\ and Cq must be independent of the plate separation L. Now 
the two limits L — > and L — > oo can both be regarded as describing a single plate alone in 
the world (as far as the outer regions are concerned). For the stress-energy tensor outside 
the plates to behave in the same way in these limits, we must take C\ — 0. The inward 
pressure on each plate is then 3C /-L 4 , so we would deduce C > from an attractive 
force. Therefore, the energy density between the plates, —Cq/L 4 , is negative, and we have 
deduced that WEC is violated. The full computation of the stress-energy tensor in QFT 
leads to the values C\ = and Co = 7r 2 /720, replicating the Casimir force formula]^] 

5 The reader might wonder how knowledge of a force, which is obtained from a difference of stress- 
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Figure 1: A spacetime plot of the energy density in a vacuum + 2-particle superposition 
state [32J. Dark areas represent negative values. 

One of the striking features of the Casimir result is the small magnitude of the leading 
coefficient 7r 2 /720 = 0.0137.... Indeed, the central message of these lectures is that the 
energy conditions are in various ways 'almost satisfied': violations are either small in 
magnitude, or short-lived, or when they are not, require disparate scales, highly noninertial 
motion, or large positive energies somewhere in the system. Indeed, it has been argued 
that the classical energy conditions might be regarded as holding in an operational sense, 
once on takes account of the positive energies present in apparatus used to produce and 
detect negative energy densities [55] . 

Nonetheless, it is clear that there is no possibility of insisting on pointwise energy 
conditions in QFT. To gain some insight into what might be possible, it is helpful to note 
that the energy densities of states formed as superpositions of the vacuum and two-particle 
states (like ip a above) tend to form fringes reminiscent of interference patterns. An example 
is given in Fig. [T] from which it can be seen that the fringes are spacelike in character; any 
timelike observer meets alternating positive and negative values and cannot 'surf a trough 
of negative energy density. This suggests seeking constraints on local averages of the energy 
density along timelike curves, and that is precisely what we will do. 

energy tensors, has apparently allowed us to compute an absolute stress-energy tensor. The answer is 
that the tracelessness of the stress-energy tensor has 'smuggled in' the extra information. The argument 
above therefore presents a choice of either accepting violations of the WEC or abandoning the standard 
electromagnetic stress-tensor based a conformally invariant Lagrangian. 
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1.3 An example of a QEI and its consequences 

The massive Klein-Gordon field in 4-dimensional Minkowski space obeys the following 
bound [22] 

i r°° 

(T 00 )Jt,x) \g(t)\ 2 dt > -— J \g(u)\ 2 u 4 Q 3 (u/m)du (1) 

for any smooth compactly supported g, and all Hadamard states^] u, where Q3 : [1, 00) — > 
IR + is defined by 

and obeys < Qa(x) < 1 with Qs(x) — > 1 as x — > 00. 



Figure 2: Q 3 (x) 

In the m = case, the bound simplifies, and actually gives a bound valid for all m > 



(Too^^^^t)^^-— / \g"{t)\*dt. (3) 



These bounds will be derived in Sec. |2.5| as a special cases of more general results. Note 
that 

• The left-hand side depends on the quantum state u, while the right-hand side is 
state- in dependent. 

• The bound is known not to be optimal. 

• The bound requires a certain degree of smoothness in g. In four-dimensions, it 
remains valid if one take g to be an element of the Sobolev space W 2,2 (1H), i.e., g, g' 
and g" are required to exist (in the distributional sense) and be square-integrable. But 



6 We will say more about these states later, but for now it is enough to know that they form a large 
class of physically reasonable states. 
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the bound does not apply to g with lower regularity, in particular, to discontinuous 
g. By 'sharp switching', one can trap arbitrarily large negative energy densities. Of 
course, no physical device is capable of instantaneous switching, as a consequence of 
the uncertainty principle. 
The QEIs contain a lot of information, as we now show. 

Scaling behaviour Put g T (t) = T^ 1 ^ 2 g(t/r). Then the bound, applied to g T is 

i r 1 r°° 

- / (T 00 )Jt,x) g{t/rfdt > -^-j / \g"(t)\ 2 dt, 



which, in the short sampling time limit r — > 0, is consistent with the fact that the expec- 
tation value of energy density at a point is unbounded below, and in the limit r — > oo 
gives 

liminf / (T 00 )Jt,x) g{t/rfdt > 
for any Hadamard state u, so the WEC holds in this averaged sense (known as AWEC). 



In Sec. 2.5 we will see how these results are modified for noninertial trajectories. 



Bounds on the duration of negative energy density Suppose that (T o) w < P for 
some interval t G [to, to + r] of time. Then, for any g G Co°((io, to + r )), 

p J \g(t)\ 2 dt > J (T 00 }Jt,x) gitfdt > -^J°° \g"(t)\ 2 dt. 

Rearranging and integrating by parts twice, this says that 

(g I 9"") . 1R 2 

— : — > — 107T O 

(g\g) 

for all such g ^ 0, where angle brackets denote the standard L 2 inner product. But the 
left-hand side can be minimized over g, to give the minimum eigenvalue of the operator 
d 4 /dt A on [to, to + t], with boundary conditions at each end corresponding to vanishing of 
the function and its first derivative]^] The upshot is that 

P > (4) 

where the numerical constant C ~ 3.17. Turning this around, in any time interval of 
duration r, the energy density must at some instant exceed — C/r 4 . Tighter results may 
be obtained for massive fields |16|. 



7 These boundary conditions emerge from some Sobolev space analysis 
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Quantum interest Developing this theme, the QEI can be regarded as asserting that 
the differential operator 

where p(t) = (Tqq) (t,x), is positive on any open interval of K, with vanishing of the 
function and first derivative at any boundaries. (A more precise formulation is to say that 
the Friedrichs extension of the above operator defined on the dense domain C£°(I) C L 2 (I) 
is positive for any interval I - the boundary conditions may be deduced from this [36J.) 
This leads to quite substantial restrictions on the possible form of p. 

For example, suppose that p has an isolated pulse, i.e., pit) = on [t\ — Ti,ti] and 
[t2,t 2 + r 2 ] with ti < t 2 and Ti,r 2 > 0. Choose a test function g e Cq°((£i — Ti,t 2 + r 2 )) 
that equals 1 on t 2 ] • Then the quantum inequality gives 

t2 p(t)dt>--^( T \g"{t)\ 2 dt + [ t2+T2 \g"(t)\ 2 dt 

and we can optimize over g to give 

* 2 A A A 



where 



167i 2 Tf 16it 2 t 2 8tt 2 min{ri, r 2 } 2 



A = inf / |<?"(t)| 2 dt 

9 JO 



with g restricted to smooth functions equal to 1 near t — and near t — 1. This amounts 
to an Euler-Lagrange equation </'" = with p(0) = 1, g'(0) = g(l) = g'iX) = 0. The 
solution g(t) = l- 3t 2 + 2t 3 gives A = 12, so 



min{ri, t 2 } 6 \ p(t) dt > 



2 



3 



2tt 



2 - 



which gives a nontrivial constraint on the extent to which a pulse (of any shape) can be 
isolated if the integral is negative. In particular, if p is compactly supported, it can only 
be compatible with the QEI restrictions if it has nonnegative integral (another version of 
AWEC). Abreu and Visser [1] have also shown that if p is the energy density compatible 
with the quantum inequalities and J p = 0, then p = 0. 

Ford and Roman [IS] first described this sort of behaviour with a financial analogy: 
nature allows you to 'borrow' negative energy density, but you must 'repay' it within a 
maximum loan term. Moreover, (excluding the case of identically zero energy density) the 
amount repaid must always exceed the amount borrowed. This is the so-called quantum 
interest effect - one may also show in various ways that the interest rate diverges if one 
delays payment towards the maximum loan term. The argument above (which is new) is 
based on [36], further developments of which can be found in [U [82] . A slightly earlier 
proof of some aspects of Ford and Roman's Quantum Interest Conjecture can be found 
in [73], but this is not as quantitative in nature as the bounds given here. 
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Application: A priori bounds on Casimir energy densities Experiments con- 
ducted in a causally convex globally hyperbolic region ought not to yield any information 
regarding the spacetime geometry outside the region. This insight has been used to anal- 
yse the Casimir effect for a long time [64] and has also been at the root of much recent 
progress in QFT in CST following the work of Brunetti, Fredenhagen and Verch [9]. Here, 
we combine it with the quantum inequalities; what follows is based on [31] (see also |20j). 




L 



Figure 3: Spacetime diagram of the Casimir plate set-up. 

Consider a region between Casimir plates at z — ±L/2 in otherwise flat spacetime and 
an inertial trajectory parallel to the plates. Let I be the distance from this trajectory to 
the nearest plate in the t — surface. Then no experiment conducted along the trajectory 
in a time interval of less than 11 can possibly know about the existence of the plates; 
it should be as if the experiment was conducted in Minkowski space. In particular the 
quantum energy inequalities apply, and [as the energy density is supposed constant along 
the trajectory] Q gives an a 'priori bound 

C 3.17 0.20 

By comparison, the known value of the Casimir energy density for the massless minimally 
coupled scalar field is 

7T 2 7T 2 3 -2cos 2 (vr^/L) 

00 ~ ~1140L 4 ~ 48L 4 " cos 4 nz/L 

which ranges between 3-7% of the bound as z varies in [— L/2,L/2\. 
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However, the a priori bound is valid even in situations where exact calculation is dif- 
ficult/impossible; it also applies to all stationary (Hadamard) states of the system. This 
partly answers the question (often emphasised by Ford): why are the Casimir energies so 
small? They are constrained by QEIs, which already gives a small leading constant in 
front of the l/£ 4 one might expect on dimensional grounds. There remains an interesting 
question as to why the Casimir energy density is a comparatively small proportion of the 
allowed bound. 



1.4 Some history and references 

The study of QEIs began with a 1978 paper of Ford [10], in which he argued that a 
beam of negative energy (described by a pure quantum state) could be used to cool a 
hot body and decrease its entropy. Ford argued from the macroscopic validity of the 
second law of thermodynamics that violations of the energy conditions must be governed 
by bounds of uncertainty principle type. That was borne out in subsequent derivations of 
quantum inequalities by Ford in a series of papers written in conjunction with Roman and 
Pfenning jJH SSj, H7J [TOj [72j |Hj, concerning Minkowski space and some static spacetimes. 
These papers established lower bounds on weighted averages of the energy densitjj^jof the 
scalar and electromagnetic quantum fields along a static trajectory, where the weight is 
given by the Lorentzian function f(t) = r/(7r(t 2 + t 2 )). Here r sets the timescale for the 
averaging. For example, the massless scalar field in 4-dimensions obeys a bound 

T(Tao(t,x)) u 3 



n(t 2 + r 2 ) - 32tt 2 t 4 

for all sufficiently nice states u and any r > 0. 

The first QEI for general weighted averages was derived by Flanagan [39] for the special 
case of massless quantum fields in two-dimensional Minkowski space. His argument forms 
the basis of a general argument for two-dimensional conformal field theories [26J that will 
be discussed in Sect. 14.11 



The bound discussed in Sect. 1.3 was derived in [22J for the scalar field of mass m > 
in Minkowski space of arbitrary dimension and for averaging along inertial curves with 
general weight functions of sufficiently rapid decay. This was generalized to some static 
spacetimes [35] for averaging along static trajectories. With some modification, the method 
also applies to the electromagnetic [69J, Dirac [27J and Rarita-Schwinger fields [88J. The 
general approach of [22J (somewhat rephrased) formed the basis for the first fully rigorous 
QEI [T7] for the scalar field, which was also much more general than the previously known 



results. We will discuss that argument in Sect. |2.4[ Generalizations to the Dirac and 
electromagnetic fields are also known [371 H21 EQ] . 

There is a significant literature on the theory and applications of QEIs — reviews can be 
found in [TgJ [T9J [75] and the recently published [T5] gives a popular but nonetheless careful 

8 In fact, the earliest papers consider negative energy-momentum fluxes, but the energy density soon 
became the main object of interest. 
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account. I mention only two aspects here. First, QEIs place significant constraints on the 
ability of quantum fields to support wormholes or other exotic spacetimes, if the fields are 
assumed to obey a QEI similar to those found for the free scalar fields [7H HHl [33] . Second, 
the link between QEIs and thermodynamics, which originally motivated Ford [3D], can be 
pursued abstractly (in a setting that includes the scalar field) [38] . 

2 Some methods of Quantum Field Theory in Curved 
spacetime 



The QEI studied in Sec. 1.3 can be derived directly by fairly elementary means [22]. 
However, it is also a special case of a rather general QEI, whose proof will be our main 
goal. To achieve this we introduce the algebraic formulation of QFT in CST; the completion 
of the proof will also need the tools of microlocal analysis in a subsequent lecture. In terms 
of literature, [2] contains much relevant material, while [87] emphasises a slightly different 
version of the algebraic approach (and does not cover microlocal analytic methods). Some 

(although emphases differ) which contains a more broadly 



of this material is based on [21 
based account of QFT in CST.^1 

Throughout, let M be a globally hyperbolic spacetime, understood to comprise a 
(smooth etc) n-dimensional manifold with time-oriented Lorentz metric and such that 

• there are no closed causal curves 

• for any points p, q, the intersection of causal futures/pasts J + (p)r\J~(q) is compact 
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2.1 The Klein-Gordon field 

We will study the formulation of the real scalar field, defined by Lagrangian density 

= \p 9 (^(V a 0)(V b 0) - (m 2 + ^)0 2 ) , 

where p g is the density induced by the metric g a b of M, R is the Ricci scalar and £ is 
a dimensionless coupling constant. The case £ = is known as minimal coupling and 
£ 7^ as non-minimal coupling. In the special case m = 0, £ = (n — 2)/(4n — 4), the 
action exhibits conformal invariance, because the Lagrangian density is unchanged under 
the simultaneous replacements 

9ab -)• g ab = tt 2 g ab 0^0 = tt 1 -^ 2 ^ 

for any smooth positive function Q, i.e., = «£g [</>]• This value of £ is accordingly 

called conformal coupling. 

9 I am aware of a number of misprints and minor errors in |21) . which I hope to correct in due course. 
10 Thc definition given in [531 US] also requires strong causality, but this is a consequence of the other 
two conditions and can be dropped 4, Thm 3.2]. 
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This field equation derived from this action is the Klein-Gordon equation 
P0 := (D g + m 2 + £R)<f> = 0, where D g = g ab V a V b} 
and the stress-energy tensor, obtained by varying the action with respect to the metric, is 

Tab = (V a y?)(V fc y?) - )^g ab g cd (y c v)(y d y) + -m 2 g ab (p 2 
+ f {g ab D g - V a V fe - G ab ) 2 , 

where G ab is the Einstein tensor. Note that the effect of the coupling constant can be seen 
in the stress-energy tensor even where the metric is Ricci flat, even though the ^R<p 2 term 
in the Klein-Gordon equation vanishes in such situations. 

The Klein-Gordon field is well-posed on globally hyperbolic spacetimes, for which we 
refer to the thorough and clear presentation of [3]. For our purposes, the main result is: 

Theorem 2.1 If M is globally hyperbolic then, to each f G C^°(M) there exists ± G 
C°°(M), with supp^ C J ± (supp f), solving the inhomogeneous problem 

= f, (5) 

Moreover, <p + /~ is the unique (distributional) solution to (|5| whose support is past/future- 
compact (i.e., the support has compact intersection with every set of the form J T (p)). The 
maps 

E ± : C °°(M) — ► C°°(M) 

are linear continuous mappings, where C^°(M) and C°°(M) are given their standard 
topologies. 

Due to the support properties, E~ (resp., E + ) is called the advanced (resp., retarded) 
fundamental solution (or Green function). In the special case where / = Pf for some 
/' G C^(M), we note that /' is both past and future compact, so /' = E ± f by uniqueness. 
Hence we have 

E ± Pf = f 

together with the initial property PE ± f = f. 

The advanced-minus-retarded fundamental solution E is defined by E = E~ — 
E + . (Warning: some authors use retarded-minus-advanced, or label retarded and advanced 
the other way round! Furthermore, in the — h ++ signature, the fundamental solutions 
to (D g — m 2 )(f) = f are minus the fundamental solutions we use; e.g., Wald's A jHZ] is our 
—E~.) Clearly <p = Ef is a smooth solution to the homogeneous equation P(j) = 0, but 
we also have an important result (cf. [3J, Thm 3.4.7]) that summarises a number of key 
properties in a compact form. 
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Theorem 2.2 The following is an exact sequence (that is, the image of each map is pre- 
cisely equal to the kernel of the next): 

_> c?{M) ^ C™{M) ^ Cf c {M) A Cf c {M) (6) 

where C^ C (M) denotes those functions in C°°(M) with support contained in J(K) = 
J + (K) U J~(K) for some compact K . 

Remark: The support of any function F G C^ C (M) has compact intersection with any 
Cauchy surface. But it is not the case that a smooth function whose support has compact 
intersection with each leaf of a given foliation of M by Cauchy surfaces is necessarily in 
Cgg(iVi")P] Unfortunately the literature contains many references to functions 'compactly 
supported on Cauchy surfaces' that would be more accurately rendered as 'in C^ c (My. 
Proof: The equalities EPf = E~Pf- E+Pf = f - / = and PEf = PE~f - PE+f = 
f — f — for / G C£°(M) are immediate, so each image is certainly contained in the 
kernel of the following map. For the reverse inclusions, we observe that 

• if Pf = with / G C^°(M) then / = E + = by uniqueness of past-compact 
solutions; 

• if Ef — with / G C^(M) then E + f = E'f, which shows that E + f is supported 
in the compact set J + (supp/) fl J~(supp/) and hence / = PE + f G PC^(M); 

• if P<f) = with G C^ C (M) we argue as follows. Choose any two Cauchy surfaces 
H ± , with S + C and a smooth function \ with \ — 1 in J~(X~) and x = 
in J + (S+). Then 

f = Px<P 

is compactly supported (in fl As x4> has future-compact support, 

X<t> = E~f. But (x~ 1)0 has past-compact support, and P(%— 1)0 = /, so (%— 1)0 = 

Subtracting, (f> = E~f- E+f = Ef. 
Note: This shows that if O is any open neighbourhood of a Cauchy surface, then any 
solution may be expressed as Ef for some / G C£°(0). 

□ 

2.2 Phase space 

The symplectic space Our phase space consists of all real-valued solutions with SC 
support 

SoIr(M) := {0 G Cf c {M- E) : P0 = 0}. 

11 For example, in four dimensional Minkowski space, the set U^ 1 {l/n} x B n , where B n is a ball of 
unit radius, centred at (An, 0, 0) € K 3 has compact intersection with each t — const hypersurface, but is 
not contained in J(K) for any compact K. 
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However, it is rather convenient to work with its complexification, i.e., the space of complex- 
valued solutions 

Sol(M) := {0 G Cf c {M) : P0 = 0}. 

In view of the exact sequence (|6| and the first isomorphism theorem for vector spaces, this 
may be reformulated as 

Sol(M) = EC™(M) S C^{M)/PC^(M). 

Let us write E for the isomorphism E : C™(M) / 'PC^(M) — > Sol(Af), which has action 

E(f + PC™(M))=Ef. 

The pairing C°°(M) x C °°(M) -> C, 

(0,/)^ / 0/dvol 

JM 

clearly induces a pairing Sol(M) x (C CO (M)/PC CO (M)) -> C by 

(<M/]> = / 0/^vol 
Jm 

(we use the formal self-adjointness of P here). We may now define a bilinear map on 
Sol(M) by 

c^i) 02) = (02, -S~ 1 0i), (note reversal of order!) 
which evidently has the properties that 

<r(Ef, 0) = (0, E^Ef) = (0, [/]) = / 0/dvol (7) 

and hence 

a(Ef\,Ef 2 )= [ (£/ 2 )/idvol^ £(/!,/.,). 
Jm 

An easy calculation shows that 

CT(01, 2 ) = y (0lV n 2 - 2 V Tl 1 ) rfS (8) 

for any smooth spacelike Cauchy surface S with future-pointing unit normal vector n, 
from which it is clear that a is antisymmetric and, moreover, is the standard symplectic 
form for the Klein-Gordon system. To prove ([§]), write 0i = Ef for some / G Cq°(M) 
supported to the past of S. Given the definition of E and the support properties of E ± , 
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h = (E- - E + )f = -E+f on S, so 

RHS of (§ = y (-E + fV n (/) 2 + <p2^ n E + f) 



/ V a {<j> 2 V a E + f - E+fVafo) dvo\ 
Ji-(s) 

[ {foPE+f - (E + f)P<f) 2 ) dvo\ 
Ji-(s) 



/-(E) 

cj) 2 fdvo\ = a(0i,0 2 ) 



M 



using the divergence theorem in conjunction with the fact that E + f has past-compact 
support, and supp/ C together with Eq. 

The map a is evidently weakly nondegenerate, in the sense that if cr(0', 0) = for 
all 0' G Sol(M), then, putting 0' = £?/, we find / 0/dvol = for all / G C °°(M) and 
hence = 0. 

As mentioned, a is the standard symplectic form for the Klein-Gordon field. To make 
our conventions more explicit, we observe that the covariant momentum conjugate IP to 
the field is defined by the functional derivative^] 



Pg 0\/ a <P 

SO 

^(0i,0 2 ) = / (0il^-0 2 n?)n a d£. 

JT, 

(This corresponds to the convention that the symplectic form, in finite dimensions, may 
be written a = dq l A dpi in terms of canonical coordinates q l and momenta pj.) 

The upshot is that SoIr (Af), equipped with (the restriction of) a is a weakly nonde- 
generate symplectic space, while So\(M) equipped with a and complex conjugation, is its 
complexification. 

Classical observables and Poisson brackets Classical observables are functions on 
this phase space: for example, every / G Cg°(M;R) defines an observable Ff which acts 
on solutions G SoIr(JW) by 

F f (<j>) = / dvol g {p)<f>(p)f(p) = a{Ef, ( f>). 

J M 



2 The meaning of this expression is that 

/ n a w a dvoi= — f &(<t>,v<j> + \w) 

Jm Jv 



for every smooth compactly supported covector field w a , and any relatively compact open subset V con- 
taining suppw. 
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(By weak nondegeneracy, the last equality shows that there are enough such observables 
to distinguish elements of Sol K (M).) We may observe that the Ff depend linearly on /, 
and that some of them vanish identically: 

F Pf (<i>) = <T(EPf,<i>)=0 

for any / G C °°(M;R) and <p G Sol R (M). 

If (&,u) is a finite-dimensional (real) symplectic manifold, the Poisson bracket of two 
smooth functions F,G G C°°(^) is given in terms of the exterior derivatives of F and G 
by 

{F,G}(p) = dF(dG%, (9) 
where {dGy\ p G T p ^, which is the Hamiltonian vector field induced by G, satisfies 

tu p ((dG)%,v)=dG\ p (v) (10) 

for v G T p ^P(M) according to our convention for the symplectic formj^] In particular, if 
3? is a vector space (regarded as a manifold with T v & = and F and G are linear 
functionals on then dF\ p (v) = F(v), etc, so the Poisson bracket-a function on phase 
space-is a constant, 

{F, G} = F{G S ), where uj(G\ v) = G{v). 

Although infinite-dimensional manifolds require care, these formulae will be enough for our 
purposes. With G = Ff, we know that 

Ff(<f>) = a(Ef,<f>) 

so we may take Fj = Ef (and there is no other solution, by weak nondegeneracy). Hence 
{F h ,F fa } = F fl (Ef 2 ) = [ f 1 Ef 2 dvo\ = E(f 1 J 2 ). 

J M 

It is worth observing that our class of observables JP(M) := {Ff : / G C£°(M] R)} is, 
itself, a copy of the phase space, when equipped with the Poisson bracket as the symplectic 
form: the map Ff i-> Ef is easily seen to be a symplectic isomorphism. So this class of 
observables provides a complete description of the underlying dynamical system. 

2.3 Algebraic formulation of the quantum field theory 

The algebraic approach is actually nothing but Dirac quantization, but without requiring 
quantized observables to act on a Hilbert space in the first instance. 

13 Our convention for Poisson brackets then amounts to 

dF dG OF dG 



G} ^ ^ ^ 



in canonical coordinates in the finite dimensional case. 
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Dirac quantization Applying Dirac's quantization prescription to the classical observ- 
ables Ff, we seek (at least formally) self-adjoint operators 14 Ff (/ G Cq c (M;M.)) obeying 
the same algebraic relations as the Ff, but with the standard replacement of Poisson 
brackets by commutators 

\F f ,fy]=i{F^F r } = iE(f,f)l (11) 

(note that any constant function is quantized as an appropriate multiple of the unit 1). As 
quantizations of the classical smeared fields, the Ff are interpreted as smeared quantum 
fields. In particular, when the supports of / and /' are spacelike-separated, Ff and Ff 
should commute, reflecting the Bose statistics of a spin-0 field. 

It is also convenient to permit smearings with complex-valued functions. Accordingly, 
we define 

$(/) = F^r f +tF\^ f 

for / G Cq°(M), dropping the hats from now on and seek to implement the following 
relations): 

• f i — y 3>(jf) is complex-linear; 

• $(/)* = $(/) for all / G C °°(M) 

• $(P/) = for all / G C °°(M; E) for all / G Cg°(Af) 

• [$(/), $(/')] = z£(/, f)l for all /, /' G C °°(M). 

This may be done by invoking a unital *-algebra with abstract elements <&(/) (/ G C^°(iVf )) 
as generators, subject to the above relations. We denote it A(M). (The only risk is that 
A(M) might be trivial, but it is not: as a vector space it is isomorphic to the symmetric 
tensor vector space 

oo 

r (Sol(M))=0Sol(M) " 

n=0 

over the solution space Sol(JVf) on M. In fact, A(M) is also simple, so one could not 
impose additional relations without it collapsing to the trivial algebra.) 



States and the GNS representation Self-adjoint elements [A* = A] of A(M) should 
play the role of observables. However, this is rather empty without a rule for turning 
observables into expectation values, in other words, notion of a state. 

Definition 2.3 A state on A(M) is a linear map u : A(M) — > C obeying 

= 1 normalisation 
\/A G A(M), u(A*A) > positivity. 

14 We will write hats on top of operators only in this section. This should not be confused with the 
notation for a Fourier transform used later. 
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Expectation values 



Wn(/l, • • • , / n ) = w($(/i)$(/ 2 ) • • • $(/«)) 



are called n-point functions. It is clearly sufficient to specify the n-point functions to fix 
u>. The algebraic relations in M have implications for the n-point functions: for example, 

W 2 (/l, / 2 ) - w 2 (/ 2 , /i) = w($(/i)$(/ 2 ) - $(/ 2 WO) = (A, / 2 )1) = / 2 ), 
and 

w 2 (P/i, / 2 ) = w($(P/0^ C/0) = = waCft, P/ 2 ) 
while positivity of the state implies directly that 

u, 2 (J, /) = u, WW)) = ^ W) W) > 0. 

Thus u)2{x, x') is a bidistribution of positive type that is a bisolution to the Klein-Gordon 
equation and whose antisymmetric part is ^E(x,x'). 

Perhaps reassuringly, given a state u we may regain a Hilbert space setting using 
the GNS construction (Gel'fand, Naimark, Segal) which gives a Hilbert space J%j, a 
dense domain @ w C Jf%j, a representation 7r w of *4(M") as (generally unbounded) operators 
defined on and a distinguished vector G ^ such that 

(n u \ir u (A)n u )=«;(A) 

for all A G *4(iVf). However, we will not develop this here; see, e.g. [52]. 

Hadamard states The algebra A(M) admits rather too many states and it is necessary 
to select a 'physically reasonable' subclass. We consider Hadamard states which are 
states whose 2-point functions are distributions and take a specific form for near-coincidence 
of the points. The Hadamard class was precisely described in [66]; the definition given there 
is rather involved, but the rough idea (in four spacetime dimensions) is that when x and 
x' lie in a common causally convex geodesic normal neighbourhood, one should have 

u 2 (x, x') = U 2 {X \ X ' ] + V(x, x') log(a + /£ 2 ) + W(x, x') (12) 
Aix z a + {x, x ) 

where U, V and W are smooth, a the signed square geodesic separation of x and x', 



taken to be positive for spacelike separation and the notation f(&+) indicates a certain 



regularization of f(a) (the Minkowski space case is given explicitly below). The parameter 
i is a length scale, necessary for dimensional reasons. The functions U and V are defined 
using £, the local geometry and the Klein-Gordon operator, along with the condition that 



15 Note: Some authors, including [13] . use a for half of the signed squared separation; our convention 
follows e.g., [66] . 
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U(x,x) = 1, which allow U to be identified as the square root of the van Vleck-Morette 
determinant, 



l(i 



A(x, x) 



det[-d^d^\a{x,x')} 



^Jg{x)g(x') 

and V to be expressed as series in a. (In general, the series for V does not converge but there 
are various standard work-arounds that I will not discuss here.) All the state-dependent 
information is contained in W. There is a much cleaner definition of the Hadamard class 



in terms of microlocal analysis - see Sec. 3.2 



The motivation for Eq. ( 12 ) is that it makes the singular part of u 2 as much as possible 
like the leading behaviour of the Minkowski vacuum 2-point function, which, for mass 
m > is 

u 2 (x, x') = lim F([-(t - t' - it) 2 + \x- x'\ 2 } 1/2 ), 

e— !>0+ 

where, for m > 0, 

Tfl I 7J1 2 

F(z) = -—-K x {mz) = — - + — \og(mz) + 0(1) 

47T 2 Z 4% 2 Z 2 87T 2 

while F(z) = l/(4ir 2 z 2 ) for m = 0. Hence 



1 Tfl 2 

uj 2 (x,x') = — — + — -\ogm 2 a + + ... 



where f(cr + ) is the distributional limit 



f(a+) = lim f(-(t - t' - ie) 2 + \x - x'\ 2 ). 

The major consequence of the definition is that the difference of two Hadamard 2-point 
functions is smooth. 

Quantities like the Wick square and stress-energy tensor can be defined by normal 
ordering relative to some reference state u R using a point-splitting prescription, e.g., 

(:$ 2 :(x)) w = (u2-w?)(x,x) 

One can do without the reference state if, instead of w-f , we subtract a local Hadamard 
parametrix, i.e., an expression of the form of the RHS of (12), but with W determined by 
local geometry rather than a reference state. Actually, there are remaining freedoms in 
W that give finite renormalisation freedoms; we suppose that some choice has been made 
and denote the resulting object by & 2 en . The procedure is described in j87j with particular 



16 See Q]Jl §1], modulo change in notation. Here the derivatives are partial derivatives in some coordinate 
system, and g = \ det g^ v \ in the same coordinates. Exercise: check that this is a bi-scalar quantity. 



23 



reference to the stress-energy tensor j^] see |58j for a much more far-reaching development. 
Of course 

with normal ordering relative to uo R . 
2.4 The QEI derivation 

Let 7 : IR — )■ M be a smooth timelike curve, with proper time parameterisation. Let Q be 
any partial differential operator with smooth real coefficients. We consider the quantity 
:(Q&) 2 :, with normal ordering performed relative to a reference Hadamard state u R , and 
seek a lower bound on 

r dr\g(r)\ 2 (:(Q^:Ul(r)) 



for g £ C£°(K) and any Hadamard state u. 

To this end, we introduce a point-split quantity 

G(t,t') = (Q<S>( 1 (t))Q<S>( 1 (t'))U 

and write G R for the same quantity evaluated in the reference state. Both G and G R are 
distributions, but their difference F = G — G R is a smooth function, which is symmetric 
[as both G and G R have equal antisymmetric parts] and whose diagonal r' = r gives 

F(r,r) = (:(g$): 2 ( 7 (r))) w . 

Then for any real- valued g £ Cq°(M.) we compute 

dr\g(r)\ 2 (:(Q^:)Mr))= [ dr\g{r)\ 2 F{r,r) 



^ d ^9(r)g(T')e- ia ^F(T,r') 



(inserting a (5-function to 'unsplit' the points) 

"°° da 



F(g_ a ®g a ) 

oo ^ 

(thinking of F as a distribution and writing g a {r) = g(r)e' iaT ) 

f°° da 

= / — F{g-a®g a ) 
Jo 71 

17 As a sketch: Let T be a differential operator that maps smooth functions on M x M (or local subset 
thereof) to smooth bi-covector fields, with the property that T(4>® (j>)(x,x) is the classical stress-energy 
tensor of any Klein-Gordon solution <f>. Applying T to cj 2 — H , where H is a local Hadamard parametrix, 
and bringing the points together, we obtain a rank-2 covariant tensor field x t— > (T(w 2 — H))(x, x). It turns 
out that although this tensor field is not necessarily conserved, the problem can be fixed by subtracting a 
local geometrical term of the form Qg a b, and can be avoided altogether by a clever choice of T [57] , 
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using the symmetry of F, and hence F(g_ a g) g a ) = F(g a <g) g_ a ) to make the final step. 
As g is real- valued, we have g- a (r) = g a {j) and obtain 

dr \g(T)\ 2 (:(Q^:)Mr)) = f°° —Ffa ® g a ) 



71 

da 



>- I —G R (g^®g a ), (13) 



o 



71 



where we have used F = G — G R and the positive type property of G, which it inherits 
from oj2- The positive type property also tells us that the integrand in the final expression 
is pointwise positive in a. The result may be generalised to complex- valued g G C^°(R) 
simply by applying the above argument to the real and imaginary parts separately. 

This derivation provides a quantum inequality on :(Q$) 2 : and hence on any other 
quantity that can be expressed as a finite sum of such quantities. In particular, it applies 
to the energy density of the minimally coupled scalar field. Note that the bound depends 
only on the reference state ujr together with 7 and g. We summarise with a theorem 

Theorem 2.4 Let M be any globally hyperbolic spacetime, Q be any partial differential 
operator with smooth real coefficients, 7 be any smooth timelike curve in a proper-time 
parameterization. For normal ordering performed relative to any Hadamard reference state 



uj r , the inequality (13) holds for all Hadamard states u of the real scalar field and all 
geC?(R). 

However, there are two important questions that must be resolved to complete the 
proof of this result: 

• Is it legitimate to restrict the differentiated two-point function to the world-line, as 
we did in defining G7 



• Is the final integral in (13) finite? (If not, then the bound would not be of much 
interest.) 

The (affirmative) answers to these questions require a more in-depth understanding of 
Hadamard states than we have previously given: namely, using some tools of microlocal 



analysis, which are developed in Section 3.1| However, the reader who does not wish to 



delve into the details should at least note that neither is simply a matter of fine precision 
because 

• the first question would be answered negatively for a null trajectory and indeed there 
is no QEI bound in this case 



one may alter the derivation above slightly to yield a bound in which the final integral 
is taken over the negative half-line and diverges. 



Remarks: 
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1. We could equally take averages of other classically positive contractions of T a b along 
the timelike curve, e.g., contracted against a null vector or possibly differing future- 
pointing causal vector fields, to obtain QNEI [32J, QDEI etc. 

2. Variants exist for averages over suitable Lorentzian submanifolds, instead of timelike 
curves (see, e.g. ITT: i . 

3. One may show that no such bounds exist for smearings over spacelike surfaces (cer- 
tainly above 2-dimensions) jl3] or, as already mentioned, along null curves [32J. 

4. The argument above, and the analogous argument for the energy density, relies on 
'classical positivity' of the quantity in question. This permits a number of related 
bounds to be proven by similar methods, e.g., see [30] for spin-1 fields. Nonetheless, 
there are also QEIs for the free Dirac field [371 021 ED] despite the fact that the 
'classical' Dirac energy density is symmetrical about zero and unbounded from below. 
It turns out that the analogue of the Hadamard condition also functions as a local 
version of the Dirac sea, and restores positivity [modulo a finite QEI lower bound] 
as well as renormalising the energy density. 

5. See Sect. [5] for discussion of nonminimally coupled scalar fields and the case of inter- 
acting QFT. 



Dependence on the reference state We can rewrite the inequality (13) as 



" da 



dr\g(r)\ 2 ((Q^JMr))> j dr \g(r)\ 2 ((Q^Y T JMl(r)) - J ^-G"fc®g a ) 

using a renormalized square, rather than Wick ordering. Now, slightly heuristically, 
((Q^)rcn)ui R (l( T )) is the diagonal of a function F ren (r, r') = G R — G Ten , where G rcn is 
formed from the Hadamard parametrix (i.e., local geometry) and the operator Q. So the 
dependence on the reference state actually cancels, and we obtain 

dT\g(r)\ 2 ((Q^)i n )Mr)) > - f° — G M (fc0 Ja ). 



o 



7T 



Making this precise and quantitative takes a bit of work [33] . 



2.5 Computations in n = 4 Minkowski space for minimal coupling 

Inertial trajectory Take Q to be a partial differential operator with constant real 
coefficients, so that Qe lkx = p(k)e lkx for some polynomial p (which necessarily obeys 
p(—k) = p(k)) and adopt the Minkowski vacuum state as the reference, with 
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With this choice, the normal ordering is precisely the conventional normal ordering of 
Minkowski space QFT (and indeed, one would normally adjust the full renormalized quan- 
tity to coincide with this as well). We take our trajectory to be j(t) = (r, x) in standard 
inertial coordinates. Then 



(2vr) 3 2lo 



and so 



, d 3 k \p(k)\ 2 ~, . f d 3 k |p(fc)|= 

In the specific case of the energy density (of the minimally coupled field) measured 
along the inertial curve, the classical expression 



/x=0 

leads us to consider the operators Q = 2~ 1 / 2 d fl for /i = 0, 3 and Q = 2~ x l 2 m. Summing, 

da f d?k u 2 + \\k\\ 2 + m 2 v , 2 



dr\g(r)\ (:T 00 :) u (t, x) > - / — J — \g(u + a 



x da f°° dk , 2 2 



7T ^0 



so changing variables from fc to u = y/k 2 — m 2 , and then from (a, to) to (a + co,co): 

1 



471-3 



in 



du\g(u)\ 2 / rfww 2 ( W 2 - m 2 ) 1/2 



16?r 3 ; du \g( u )\ u Qsiu). 
where Q3 was defined in (J2|. This is the bound stated as Q. 

Uniformly accelerated trajectory Again we use the vacuum state as the reference, 
restricting to massless fields for simplicity. Here it is more convenient to work with the 
vacuum two-point function in the form 

1 

oj 2 (x,x') 



4ir 2 t7 + (x, x') 



rather than to use a Fourier representation. We consider the trajectory 

7 (r) = (eosinh(r/e ),eocosh(r/e ),0,0) 
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in inertial coordinates, where £o > is constant. This is easily seen to be a proper-time 
parameterisation of a trajectory with uniform proper acceleration ^q 1 . 

If we introduce coordinates x(t, £, y, z) = (£ sinh(r/£o), £ cos h( r /£o); y, z ) then our tra- 
jectory is £ = £o; y = z = 0. Moreover, on this curve, the vectors d/dr, d/d£, d/dy and 
d/dz form an ortho normal basis and the point-split energy density may be written as 

1 f d 2 d 2 d 2 d 2 \ 

2 \drdr' + d£d? + dydy' + dzdz' ) : ^ x ' x )- 



In terms of the above coordinates, we have 
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1 r , , , x 

u 2 (x,x') = lim — — 

e^0+ 47T 2 



(£ - £') 2 - 4^' sinh 2 ( T ^ ) + (y - y') 2 + (* - ^') 2 



and after some calculation one finds that the point-split energy density of the reference 
state restricts to the trajectory as a boundary value distribution 

3 f cy 2zc 

G r (t,t') = T(t — t'), where T(a) — lim . cosech 4 



e _ + 32tt 2 £ 4 V 2£ 

The Fourier transform may be shown to be [31] 

{) 2n \1 - e- 2 *t°" 
and similar calculations to those above give 

dr\g(T)\ 2 riH-.T ab :)Mr)) > -— / du\g(u)\ 2 T^ ,u) (14) 

1D7T J _ OQ 

for any Hadamard state u, where 

v 3 + < - 2 



T(£ ,n)=4/ cfo 



As |g(w)| 2 is even, we may replace Y(£ ,w) by \ (T(£ ,u) + T(£ ,—u)) in (14). But one 



18 A little justification is needed here, because our standard it regularisation gives 

a e = (£ - O 2 - sinh 2 (^^) + ^ sinh |- - £' sin h + e 2 + (y - y'f + (z - z'f 

in these coordinates. The important point is that a t has positive imaginary part when x and x' are null- 
separated with x to the future of x' (which implies r > r'). This is also true of the corresponding term in 
our expression for u)%, and so the alternative ie prescription is valid. 
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easily sees that 

v 3 + v^ 2 A f° ^ v s + < - 2 



T(6,,tt) + T(Co, -«) = 2T(£ , 0) + 4^ dv - 4 jT <fo ^ 



e 



2T({ , 0) + 4^" A, + (— ^ + ^^s^) 



"V - 

= 1 



4 11 

U + 



where we have combined the integrals in the second step, changing variables from u to 
— u in one of them and inserted a closed- form expression for T(£o,w). Accordingly, the 
following QEI holds for all Hadamard states u and all real- valued g G Cg°(M): 



|s(r)r 7 Y(:T a ^(7M) > / rfw 1^)1 K + 72 + 



327r»y_ DO Iav 71 V Co 2 SO^o 4 



(15) 

where «o = l/£o is the proper acceleration of the trajectory. Comparing with ([3]), we see 
that the acceleration leads to modifications to the QEI bound that are lower order in the 
number of derivatives applied to g. 



In particular, the scaling behaviour discussed in Section 1.3 is modified; we have 



1 W*>.(7<*))| 9 <*/r)l»*> " J "" 2 ll9l|2a ° 11|ljl|2Q ° 



T ' ' ^ ««/"V'W/l*w /I - 



where the norm is that of L 2 (IR). For oqt <C 1 the previous result is recovered to good 
approximation; however, for a^r ^> 1 it is the last term that dominates and, indeed, 
the AWEC fails - multiplied by r, the right-hand side diverges to —00 as r — > 00. By 
subtracting this troublesome term we can deduce that 

liminf / (j«i b (T ab )Mt)) + ^) W/r^dt > 

holds for any Hadamard state u. It is a remarkable fact that the constant negative con- 
tribution is precisely equal to the energy density of the Rindler vacuum state for the 
right-wedge {(i, x, y, z) : x > \t\} of Minkowski space. (Even though the Rindler vacuum 
does not extend to a Hadamard state on the whole of Minkowski space, it is Hadamard on 
the interior of the wedge, which completely contains the accelerated trajectory. Arguments 
based on local covariance [5] show that the Minkowski QEI along that trajectory must be 
respected by the energy density of the Rindler vacuum - see [31] for discussion and other 
examples - but it is nonetheless surprising that the Rindler vacuum saturates the QEI in 
this way.) 
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This example might suggest that a good way of 'mining' negative energy density is 
simply to follow a uniformly accelerated trajectory, when the field is prepared in (an 
approximation to) the Rindler vacuum. It is worth noting that the work required to 
maintain this motion grows exponentially with the proper time, and therefore the 'cost' in 
work done is growing much more rapidly than the 'benefit' of negative energy 'seen'. This 
seems to fit a broader pattern of adverse cost-benefit analyses in other situations where 
sustained negative energy densities may be created. 



3 Microlocal analysis and Hadamard states 
3.1 The wavefront set 

Fourier analysis provides a fundamental duality between smoothness and decay: smooth 
functions have rapidly decaying Fourier transforms, and vice versa. The fundamental 
idea underlying microlocal analysis is that decay properties of the Fourier transform of a 
distribution can be used to obtain detailed information about its singular structure. A 
general reference for this section is [61J, particularly chapter 8. 

The wave-front set Recall that our convention for the Fourier transform of functions 
is 

jn„ Ak-x 



f(k) = J d n xe lk x f{x). 

The Fourier transform of a compactly supported distribution T is, similarly, T{k) = T(efc), 
where Ck{x) = e . The Fourier transform of Schwartz distributions can be defined using 
duality, because the Fourier transform is an isomorphism of the Schwartz space to itself, 
and hence dually of the Schwartz distributions; general distributions in S 1 ' do not have 
Fourier transforms. 

The duality between smoothness and decay mentioned above is illustrated by the fol- 
lowing examples. 

a) If / e C™(R n ) then 

(l + |fc| 2m ) f(k) = |(1 + (-A) m f) A (k)\ < J <Px\(l + (-A) m f)(x)\ < oo. 

So for each N, there exists a constant Cn such that 

rr, 

/(*) 



-TTW 



(this is what we mean by 'rapid decay'.) 

b) The (^-distribution at the origin has Fourier transform S(k) = 1, which exhibits no 
decay at oo. 
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c) The distribution T G 3>'(R) defined by 

T(f) = lim f -tt$-ds 

e~^0+ J S — is 

has Fourier transform 

-Aks 



ds = 2mQ(k), 

s-ie 



which decays as k — > — oo but not as k — > +00. 

The wavefront set localises information of this type both in x-space and on the "sphere at 
00" in /c-space. 

Definition 3.1 (A) If u G @'(R n ), a pair (x, k) G R n x ((M n )*\{0}) is a regular direction 

for u if there exist 

1) G C °°(R n ) with (f>(x) ^ 

ii) a conic neighbourhood V of k in (R n )* 

Hi) constants C^, N G N 
so that 

< ITW ykeV ' NeN 



4>u{k) 



i.e., (pu decays rapidly as k — )■ 00 inV. 
(B) The wavefront set of u is defined to be 

WF(u) = {(x, k) G R n x ((]R n )*\{0}) : (x, k) is not a regular direction for u}. 

Examples 

a) If / G C°°(R n ), then WF (/) = 0. 

b) WF(£) = {(0,k) G R 2 : k ^ 0}. (Note that 05(k) = 0(0), so (x, k) is a regular 
direction for x 7^ as we may then choose with <f>(x) 7^ 0, 0(0) = 0). 

c) WF (T) = {(0, fc)Gl 2 :fc>0} (exercise!). 

The wavefront set has many natural and useful properties. For our purposes, the most 
important are the following: 

• WF(«) = u e C°°(R n ). 

• WF(Aw + fiv) C WF(u) U WF(u) for u, v G $'(R n ), A, fi G C. 
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• If P is any partial differential operator with smooth coefficients, then 

WF(Pw) C WF(ti) C WF(Pu) U CharP, 

for any u G ^'(W 1 ), where CharP is the characteristic set of P. To define the 
characteristic set, let m be the order of P, i.e., the least m G No so that P may be 
written in the form P = ^2,\ a \ <m a a (x)(iD) a where a is a multi-index. The principal 
symbol of P is the smooth function on IR n x (W 1 )* given by 

p m (x,k) = ^ a a (x)k a 

\a\=m 

and the characteristic set is 

CharP = {(x, k) £l"x (W l * \ {0}) : p m (x, k) = 0}. 

• Propagation of Singularities: WF(«) \ WF(Pw) is invariant under the Hamilto- 
nian flow generated by p m . 

• Under coordinate changes, WF and Char transform as subsets of the cotangent bun- 
dle: given a diffeomorphism ip, define (p*u by ((p*u)(f) = u(f o <p). Then 

WF(«) = {(x,£D<p\ x ) : (<p(x),Z) G WF(yAz)}; 

similarly, setting (P^f) op — P(f o <£>), we have 

(P<p)m(<p(x),g) =p m (x,€D<p\ x ). 

Here £D(p\ x is the composition of £ and D(p\ x as linear maps, i.e., the action of the dual 
map to Dtp\ x on £.) In particular, we may extend the wavefront set and characteristic 
set to distributions and partial differential operators defined on manifolds; both are 
subsets of the cotangent bundle. 

Examples: 

1. Let P be the Klein-Gordon operator P = O g + m 2 + £P on a spacetime M . The 
principal symbol is easily seen to be 

P2(x,0 = -9 ab (x)Uo 

and so the characteristic set is 

CharP = Af := {(x,£) G T*M : £ a non-zero null at p}, 

the bundle of nonzero null covectors on M. Hence the wavefront set of any (distributional) 
solution to Pu = obeys 

WF(m) C M; 
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moreover, WF(m) is invariant under the Hamiltonian evolution A h- >■ (x(A),£(A)) G T*M 
given by the 'Hamiltonian' P2(x,£) = —g ab (%)£a£b- The solution curve A i-> (x(A),£(A)) is 
such that x(A) is a geodesic [which is easily seen by noting that the 'Lagrangian' underlying 
p 2 is — \g a bX a x b } to which £(A)" is tangent and along which £(A) is parallel-transported. 

Recalling that WF (u) C Af, we may deduce that if (x, £) G WF («), then £ is null, and 
further, the wavefront set contains every point (x(A),£(A)) for A G R, where x(A) is the 
null geodesic through x with tangent and £(A) is the parallel transport of £ along x(A). 
2. Now consider Klein-Gordon bisolutions, i.e., F G 3)'{M x M) such that 

[P®\)F = (1®P)P = 0. 

Now the operator P <g) 1 has principal symbol 

p(x J e;x / ,0 = -^)e«6 

and characteristic set 

Char (P <g) 1) = (A/" x T*M) \ Z 

where Ao is the bundle of (possibly zero) null covectors on M (i.e., A/", with the zero 
covector added at each point) and Z is the zero section of T*(M x M). Similarly, 1 £g> P 
has principal symbol 

p'(x^;x',o = -g ab (x')Ca 

and characteristic set 

Char (1 ® P) = (T*M xjV )\Z 
Any bisolution F therefore has wavefront set with upper bound 

WF (F) C Char (1 ® P) n Char (P <g> 1) 

C ((A/" x T*M) n (T*M x A/" ))) \ 2 
C (A/" x A/" ) \ 2 

Pull-backs Suppose X and K are smooth manifolds and <p : K — >■ X is smooth. Given 
it G f^'(Jf), Theorem 2.5.11' in [60] constructs the pull-back <yj*ii as a distribution on F 
provided WF (u) fl = 0, where 

N v = {(<p{y), G T*X | ££V(y) = 0} (16) 

defines the set of normals of the map (p. The wave front set of the pull-back is constrained 
by 

WF {<p*u) C </WF (it) = {{y, ZD<p{y)) \ (<p{y),() G WF (it) } . (17) 

If u is smooth, the pull-back reduces to ordinary composition (p*u(y) = u(<p(y)). 
Example Let F G &{M. x Af) be a Klein-Gordon bisolution, and let 7 : R — >■ M be a 
smooth timelike curve. We wish to consider the pull-back (7 x j)*F. 
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To see that this is well-defined, we first compute the set of normals to 7 x 7 : IR 2 — y 
M x M, noting that 

(k k>) L>( 7 x 7 )| (TiT0 Q = ( 7 (r) • fc)f + ( 7 (r') • fc')f 

and therefore vanishes for all t, t' iff 7(7-) • k = j{r') • fc' = 0. Thus 

AT 7X7 = {( 7 (r), k; 7(7-'), fc') G T*M x T*M | 7 (r) • k = 7(1-') • fc' = 0} . 

Now the covectors arising in WF (F) are always null and at least one of them must be 
nonzero; moreover, no nonzero null covector can have vanishing contraction with a timelike 
vector. Thus 

WF (F) n iV 7X7 = 
and the pull-back is well-defined, with wave-front set 

WF (G) C (7 x 7 )*WF (F) 

= {(r, 7 (r) • k; r, 7(7-') • k') G T*R x T*R \ ( 7 (r), k; 7(r'), fc') G WF (F)} . 

The same is true for any distribution QF, where Q is a partial differential operator on 
M x M with smooth coefficients, because the wave-front set cannot expand under the 
action of Q. 

There are similar wave-front set conditions under which products of distributions can 
be defined. 



3.2 Microlocal formulation of the Hadamard condition 

Let us compute the wave-front set of the Minkowski vacuum 2-point function 

Consider a localising function of the form 0(27, x 2 ) = 4>i (27)02(^2) where 0j G Cq'(M). 
Then 

toW) = J j^f^M? - k)M* + k ) 

with future pointing, on-shell k. As the functions fa are smooth, their Fourier transforms 
decay rapidly as their arguments become large. The main contribution to the integral 
therefore arises from regions of k where t — k and t' + k are simultaneously small, i.e., £ 
must be near to the future pointing on-shell covector k, and £' must be near —k. Arguing 
in this way, it is not hard to see that there are open conic neighbourhoods of x IR 4 and 
M 4 x M 4 in which the integral will tend rapidly to zero as {1,1') — > 00, where IR 4 , is the 
half-space in which ±k > 0. 
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Figure 4: If t — > oo in the shaded region then (fiu 2 (l, I') — > rapidly, regardless of £'. 

Thus (xi, ki] x 2 , k 2 ) is a regular direction if either (i) (ki) < or (ii) (k 2 )o > 0. Putting 
this together with the upper bound 

WF (co 2 ) C M x J\f 

we conclude that 

WF (u 2 ) C N + x Af- , (18) 

where 

J\f ± = {(p,f) G Af : ^ is future(+)/past(-) directed}. 

Now any Hadamard state of the Minkowski theory must have the same 2-point wave-front 
set, because 2-point functions of Hadamard states differ by smooth functions. We now 
elevate this to a general principle in curved spacetimes. 

Definition 3.2 A state u> obeys the Microlocal Spectrum Condition (fiSCj^if 

WF(oo 2 ) C Af + x Af~ . 

In particular, this asserts that the 'singular behaviour' of the two-point function is 
positive-frequency in the first slot and negative- frequency in the second. We have already 
argued that the Minkowski vacuum obeys the /iSC; it is also true that ground and thermal 
states on various classes of stationary spacetime satisfy the /xSC [621 EE El]- (In relation 
to thermal states, the key point is that negative frequency contributions to the first slot 

19 The term microlocal spectrum condition was introduced in [5] with an apparently stronger definition; 
sec the remarks at the end of this section. 
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of the thermal two-point functions are exponentially suppressed, rather than being absent. 
But this is enough to get the necessary decay properties.) 

A truly remarkable fact is that the /xSC is enough to completely fix the singular structure 
of u>2 (and even more: see remarks at the end of this section). Bear in mind that even if 
two distributions have the same wavefront set, their difference is not necessarily smooth 
(WF (5) = WF (25), for instance). The following result is due to Radzikowski [73] . 

Theorem 3.3 If u and uj' obey the fiSC then 

oj 2 -u 2 E C°°(M x M) 

i.e., the fiSC determines an equivalence of class of states under equality of two-point func- 
tions modulo C°° . Moreover, the fiSC is equivalent to the Hadamard condition. 

As mentioned, it is surprising that such a result can be true. The key point is that, 
while the antisymmetric parts of u 2 and u' 2 are both equal to \%E, WF(u 2 ) is not the whole 
of WF(E), which intersects both Af + x J\f~ and J\f~ x M + . Accordingly, the singularities 
in the symmetric part must precisely cancel the unwanted singular directions in WF (E) , 
which is how the microlocal spectrum condition does, after all, fix the singular structure 
of the two-point function. 



It follows from Theorem 3.3 that all two-point functions of states obeying /iSC must 
have equal wavefront sets. The universal nature of the antisymmetric part of the two-point 
function also allows us to fix the wavefront set of the two-point function as follows. 

Lemma 3.4 If u obeys the fxSC then WF(u 2 ) = WF(E) n (A" + x Af~). 

Proof: Define u 2 (x,x') = uj 2 (x',x), so WF (w 2 ) C Af~ x A/" + by the fiSC. But, using 
successively that %E = u 2 — u 2 and u 2 = Cj 2 + iE, we find 

WF (E) C WF (w 2 ) U WF (u) 2 ) C WF (u 2 ) U WF (E) 

so, using again the fact that WF (u 2 ) C Af~ x Af + 

WF (E) c WF(u 2 ) U (AT x JV+) C (A"" x A"+) U WF (E), 

and we take intersections with A/" + x A"~ to obtain the required result. □ 

The wavefront set of E is known from work of Duistermaat and Hormander on dis- 
tinguished parametrices. This permits us to give a final form of the wavefront set of a 
Hadamard 2-point function: 

WF (u 2 ) = {(p,£;p, -C) eT*(MxM)\Z: (p,£) ~ ( P ',0 and f e Af + }, (19) 

where the equivalence relation ~ is defined so that (p, £) ~ (p 1 , £') if and only if either 

• there is a null geodesic 7 connecting p and p', so that £ is parallel to 7 b at p, and £' 
is the parallel transport of £ to p' (and necessarily parallel to 7 b at p'); or, 

• p = p' and £ = 
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Eq. (19) is the form that Radzikowski stated as his 'wavefront set spectral condition'. 
Remarks: 



1. As the following quotation, taken from the 1978 paper of Fulling, Sweeny and Wald 
|50j . makes clear, the introduction of the Hadamard condition was a spur to the 
development of the algebraic approach to QFT in CST: 

All these considerations suggest that the validity of [the Hadamard condi- 
tion] be regarded as a basic criterion for a "physically reasonable" state, 
perhaps even as the definition of that phrase. This raises the possibility of 
constructing quantum states from two-point distribution solutions of the 

field equation by a procedure of the Wightman or GNS type bypassing 

the quantization of normal modes in a Fock space. 

2. We have only discussed regularity of the 2-point function. In some references, the 
term microlocal spectrum condition is defined as a condition on all n-point functions 
of the form 

wfk) cr n , 

where the T n are particular subsets of T*M xn . This condition was introduced in [5], 
where it is also shown to be satisfied by all quasifree Hadamard states. Very recently, 
Sanders has proved that this apparently more general condition is actually equivalent 
to the fiSC in the form we have stated; and, moreover, that all states obeying /iSC 
have smooth truncated n-point functions for n ^ 2 [771 178] . One may interpret this 
as saying that all Hadamard states are 'microlocally quasifree'; it also shows that the 
class of (not necessarily quasifree) Hadamard states is precisely the 'state space of 
perturbative QFT' studied by Hollands & Ruan in [57], and previously identified as 
a plausible class of interest by Kay [65] . 

3. Finally, we mention a variation on the theme. For some purposes it is sufficient only 
to require the two-point functions to agree with a Hadamard parametrix modulo cor- 
rections in some Sobolev space, rather than modulo This leads to the microlocal 
study of adiabatic states [63J. 



3.3 Application to QEIs 



There were two issues to resolve in completing the proof of the general QEI in Sect. 2.4 



The first was to establish the validity of restricting the differentiated 2-point function to 
the worldline. Effectively we want to define 

G{t,t') = {(Q®Q)u 2 )){ 1 (t), 1 {t% 

i.e., G = (7 x 7)* {{Q <g> Q)u2)). This is well-defined by the example at the end of Sec- 
tion 3.1| moreover, 

WF (G) = {(r, 7 (r) ■ k; r', 7 (r') ■ k') G T*R x T*R | ( 7 (r), k; 7 (r'), k') G WF (u 2 )} 
cRxR+xRxR", 
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because the covector in the first slot of WF (u 2 ) is future-pointing, as is while the 

covector in the second slot is past-pointing. Here we have used both the Hadamard con- 
dition and the timelike nature of the curve in an essential way. The same results apply to 
G R , of course. 

The second question concerned the convergence of 

/ — G R (g^®g a ). 
Jo TT 

Now the integrand is 

J drdT'g(T)e~ iaT g(T')e iaT 'G R (r,r') = [(g ® g)G R ] A (-a,a) 

and this decays rapidly as a — > +00 by definition of the wave-front set, and the bound 
WF (G R ) C lx R + xlx M - . Thus we have convergence of the integral and a finite bound 
- and we also see why it would have been a bad idea to arrange the final integral in terms 
of an integral over the negative half-line. 



4 Conformal field theories 
4.1 Derivation of the QEI 

Conformal quantum field theories in two-dimensions provide examples of non-free fields 
for which quantum inequality results may be derived. The basic idea was given by Flana- 
gan for massless scalar fields. It was generalised to massless Dirac fields by Vollick [S3] 
and made into a general and rigorous argument for CFTs by Fewster Sz Hollands We 
will not emphasize analytical details here, although everything can be made precise and 
rigorous. Throughout, we work in two-dimensional Minkowski space; a general reference 

is [EE]. 

Recall that the stress tensor T ab in CFT is traceless and splits into chiral components 

T 00 (t,x) = T R (t-x)+T L (t + x) 
T 01 (t,x) = T R (t-x)-T L (t + x) 

and that the left- and right-moving chiral components Tl and T R commute and obey the 
spectrum condition 

Pl,r = J T L>R (v) dv > 0. 

The important feature of CFTs we will use is that reparameterisations of null coordi- 
nates v = t — x, u = t + x are unitarily implemented, in the following sense. Under the 
correspondence v 1— > z(v) := (i — v)/(i + v), the real-line is mapped to T \ { — 1}, where 
T = S 1 is the unit circle in C. If a reparameterisation v i-> V(y) lifts to an orientation 
preserving diffeomorphism of S 1 , then there is a unitary U R (V) s.t. 

U R {V)T R {v)U R {V)- 1 = V\vfT R {V{v)) - ^-{V,v}l, 
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where Cr is the central charge (for right-movers) and 

{v,v} = -iVnvj 0,2 1 



dv 2 y/V'(v) 
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is the Schwarzian derivative. The same is true for Ti and the unitaries Ul, Ur commute 

Consider one of the stress-tensor components, say, T(v) and let / be a smooth compactly 
supported positive real-valued function. We define 

T(f) = J T(v)f(v)dv 

and aim to show that there is a lower bound on the expectation values (T(/))^. 
The idea is to define v h- > V(v) by V'(v) = l/f(v) and set if)' = U(V)ip. Then 

f(v)(T(v)U = f(v)(U(V)T(v)U(V)-% 

= V'(v)(T(V(v)))v--^{V, v}f(v) 

so 

(T(v))*f(v)dv = J(T(V))^dV-^J{V,v}f(v)dv 
= { p h'-^ j{V,v}f(v)dv. 
Using the spectrum condition, and rearranging, we find 

(T(v)Uf(v)dv>-^ J (^v/JM) 2 dv 
for all 'reasonable' ip. 

The only problem is that the map v t— > V(v) does not lift to a diffeomorphism of S 1 . 
The resolution is to replace / by 

fneiv) =f(y)+e- (ne(v) 

with ( ne chosen so that 

• V'(v) = 1/ f ne (v) does lift to a diffeomorphism of S 1 for all n G N, e > 0; 

• J (T(v))^ ne (v) dv — > as n — > oo for each fixed ip, e; 

2n More generally, the theory contains commuting 'left' and 'right' unitary multiplier representations of 
the universal covering group of the orientation-preserving diffeomorphisms of S 1 , obeying 

U{p)U{p') = e lcB ^U{ PP ') 

where B is the Bott cocycle and c is the central charge. 
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d r^-V . r/d 



Then, for each fixed ^ 



(T(v))^f ne (v)dv > J (J^Vfne( v ) ) dv 



so taking n — >■ oo... 



2 



1Z7T n-5>oo I \ (IV 



...and e — ?• + , we obtain the desired bound 



d 



(T{v)Uf(v)dv>~J f^V/OO) *>■ 

The fully rigorous argument for this is given in [26], where an axiomatic setting is 
adopted in which all the above manipulations are justified and the class of 'reasonable' ip 
is specified. The axioms are shown to hold for CFTs constructed from unitary, positive 
energy Virasoro representations. We also proved that the bound is sharp if the theory has 
a conformally invariant vacuum. Any nonnegative / G J^(R) can be used for smearing. 

This argument is notable, partly as the first examples of QEIs for non-free fields, but 
also because it does not depend on a 'sum of squares' form of the energy density. It is also 
model-independent, applying to all unitary positive energy CFTs in one go. 

4.2 Probability distributions 

Everything said so far concerns the expectation value of the smeared stress-energy tensor or 
other similar quantities. Here, we discuss what information can be gleaned concerning the 
underlying probability distribution of individual measurements of such quantities. Again, 
CFTs provide a framework in which this can be studied for a whole class of models. The 
argument given here is taken from [23] and approaches the probability distribution through 
its moment generating function 



n=0 



IV. 



Our notation is 



G n (u n , . . . , Ul ) = (Q | T{u n ) ■ ■ ■ r(ui)fi) (Go = 1) 

g n [f} = G n (f 1 ...j) = (n\T(frn) 
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and we assume that test functions are real-valued and rapidly decaying at infinity. The 
main tool used in the argument is the CFT Ward identity [5H p. 21 



21 



n-1 



G n (u n , ...,u 1 )=} j 

3=1 



c G n - 2 (u n -i, ...,Uj,... ,ui) 



(«« - uj - toy 



djGn-^Un-i, ■ ■ ■ G n „iQ n _i, ...,U 1 ) 

2n(u n — Uj — zO) n{u n — Uj — iO) 2 



where the hat denotes an omitted variable. 

Since Go = 1 and Gi(ui) = (T(ui)) = 0, it follows immediately that 

G 2 {u 2 ,ui) 



and if we smear the Ward identity against n copies of /, we find 



n-l 



g n [f] = (n-i)g 2 [f]Qn-2[f] + J2 I n^ 

3=1 



where 



1 

*2tt 



du n duj f(u n )f(uj 
2 



u n — Uj — iO 



9jG n -i(f 



, . . . ,Uj, 



f) 



(u n - Uj - iO) 

3 

after integration by parts. Here, / * f is 



G n -l(f,---,Uj,...,f) 



du n f(u 
du 
du 
du 
du 



f( U 3 



/(«*: 



* 2n(u n — Uj — iO) J ir(u n — Uj — iO)' 

f(u n )f(Uj) f(u n )f(Uj) 



n \2ii(u n — Uj — iO) 2ir{u n — Uj — iO) 2 

f(u n )f(u 



2ii{u n — Uj — iO) 
f(u n )f'( Uj ) - f\u n )f{u 3 ) 
27t(m„ — Uj — iO) 

f{u n )f{ Ul ) - f'(Un)f(u S ) 

2n(u n - Uj) 



' + f(u n )f(uj)d n — — 

27i(u n — Uj — iO) 



(20) 



(21) 
(22) 



21 Beware, however, a misprint in Eq. (3.12a) of [5T] [©( ' should be 0( + )]. Fortunately the result given 
before (3.15) of [5T] is correct. 
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after using the Leibniz rule, a further integration by parts in one term and observing that 
the the numerator in the penultimate integrand vanishes as 0(u n — Uj) as u n — > Uj. Note 
that no boundary terms arise when integrating by parts provided / is compactly supported, 
for instance. 

To solve the recurrence relation, consider a 1-parameter family of test functions (f\(u)) 
solving 

-jY — fx * fx fo — f- 
Then the recurrence relation becomes 

g n [h] = {n- i)g 2 {h\g n - 2 {h] + ^Gn-i[fx] 

and gives a p.d.e. 

— M[fMf x ] = fj,Q 2 [fx]M[p,f x ) + Q^Mivfx] 

using C?o[/a] — 1) £i[/a] = 0. Solving, using the fact that M[/i/ A ]| M=0 = 1 and setting 
A = 0, 

M\ M f]=exp( TdX (/x-A)&[/ A ] 



a result first obtained by Haba [53J. 

In general, it is not easy to take this further. However, if / is Gaussian, we may proceed 
to a closed form result. Let f(u) = e~ u ' T j{r\pK) and make an ansatz f\(u) = A(X)f(u). 
Then 

U»*A)(«) = ^^ = ^A(«), 

so the flow equation for fx reduces to 

A'(X) = 

T 71 

The unique solution with A(0) = 1 is 

TIT 2 
7TT Z — A 

Thus Q2 [fx] = A(\) 2 G2[f] and we calculate 

c/24 

M\pf] 



1 - ^/(ttt 2 ) 

The probability distribution itself is then obtained essentially by inverse Laplace trans- 
formation: we seek P(uS) such that 



/oo 
-oo 
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and the solution is 



T(a) 



(a shifted Gamma distribution) with parameters 



a 



(3 = 7TT 2 



24ttt 2 ' 24 
which has an integrable singularity at lower limit for c < 24 




Figure 5: The probability density P(x) plotted for c = 1, where x 



TIT Ul 



One should mention that the Hamburger moment theorem guarantees that this is the 
only possible solution: noting that it is a solution, we may use it to read off the moments 
a n of P{oS) and note that they obey a bound |(?n[/]| < BC n n\ for some constants B, C, 
thereby satisfying the hypotheses of the Hamburger uniqueness theorem [79] . 

The probability distribution is clearly highly skewed. We see that the lower bound of 
the support coincides precisely with the sharp lower bound on the expectation values of 
T(f) for the Gaussian /, i.e., —c/(2Atit 2 ) - as it should for general reasons [23]. Thus the 
QEI bound, which was originally derived as a constraint on the expectation value of the 
smeared energy density, also turns out to be a constraint on the minimum value that can 
be achieved in an individual measurement of this quantity. 

The probability of obtaining a negative value is given in terms of incomplete T-functions: 



Probfw < 0) 



r(c/24,c/24) 
r(c/24) ' 



For c = 1, this results in a value 0.89 - an overwhelming likelihood of obtaining a negative 
value from a measurement in the vacuum state. In the limit c — > oo the probability tends to 
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1/2, to be expected from the central limit theorem. (Note: the above computation refers 
to just one of the chiral components of T a i, and for averaging along the corresponding 
light-ray. For Gaussian averages of the energy density along an inertial curve and c = 1, 
the probability of obtaining a negative value is 0.84.) 

It is somewhat ironic that negative energy densities, which are suppressed (in all phys- 
ically reasonable states) by the uncertainty principle expressed in the QEIs, turn out to 
occur with such high probability, in individual measurements made in the vacuum state. It 
is not known with what probability negative energy densities occur in any other state, or 
for test functions other than a Gaussian, but it would be of interest to extend the analysis 
further. Investigations of similar phenomena in four dimensions can be found in [21]. An 
application to two-dimensional dilaton quantum gravity was made in [TT] . to argue that 
the positive energy tail causes strong focussing of light cones near the Planck scale. 

5 Other directions 
5.1 Nonminimal coupling 

As mentioned above, the classical minimally coupled scalar field obeys the WEC by virtue 
of a decomposition of the energy density as a sum of squares. This is not true for the 
nonminimally coupled field, and indeed the energy density can be made arbitrarily negative 
at any given point (see, e.g., |1H] for a discussion). It turns out that this behaviour is, 
nonetheless, constrained by locally averaged energy conditions reminiscent of the QEIs [28] . 
For example, if 7 is a complete causal geodesic with affine parameter A in a spacetime M, 
and the coupling constant is £ G [0, 1/4] then there is a bound 

^dA T ab Yi b \g\ 2 > -2£^dA ||^| 2 + l -R ah i a i b \g\ 2 - Q -i) Ri 2 \g\ 2 ) <P 2 

for any solution to the nonminimally coupled Klein-Gordon equation, with corresponding 
stress-energy tensor T a b and any g G Cq°(R). In particular this result includes the case of 
conformal coupling. Note that the bound involves the field, but not its derivatives, while 
the quantity to be bounded involves field and derivatives, including some of second order. 
This inequality therefore exhibits the 'gain in derivatives' phenomenon that occurs in the 
Garding inequalities of pseudodifferential operator theory. 

The corresponding quantum theory was discussed in [29] for the case of Minkowski 
space. It was found that the nonminimally coupled field can sustain large negative energy 
densities for long periods of time. The argument is the following: the failure of classical 
WEC allows the existence of one-particle states with negative energy density near the 
origin, say. By a scaling argument these can be taken to have any desired spacetime extent, 
although the magnitude is correspondingly reduced. However, we may tensor together as 
many of these one-particle states as we wish, with respect to which the energy density is 
additive (as this is a free theory). Thus states of arbitrarily negative energy density can 
be sustained over arbitrarily large spacetime volumes. 
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However, there is a cost. The overall energy of these states is positive, and grows more 
rapidly than the scales characterising the negative energy density effect produced. So the 
production of negative energy density is inefficient in this sense. By modifying the QEI 
arguments discussed in these notes, one can establish QEI bounds for the nonminimally 
coupled field that are state-dependent [29]. In these bounds the averaged energy density 
is bounded below by state-independent terms together with terms that involve averages 
of the expectation value of the Wick square of the field in the state of interest - again 
demonstrating a 'gain in derivatives' phenomenon. If one estimates these terms using the 
Hamiltonian operator, it again emerges that the production of sustained negative energy 
densities only occurs when disproportionately large positive energies are available. This is 
still a comparatively new development; further work, it is hoped, will clarify these issues. 

5.2 Interacting fields 

There is now a fairly complete theory of quantum energy inequalities for free fields in 
globally hyperbolic spacetimes of any dimension (although optimal bounds are lacking in 
general). As we have seen, similar results hold in a large class of conformal field theories 
in two dimensions. The situation for interacting fields is more complicated, of course, and 
not so much is known. The following remarks summarise the state of knowledge: 

• One cannot expect state-independent QEIs to hold; as mentioned, these can even fail 
in the nonminimally coupled theory. Moreover, on physical grounds, we can expect 
that long-lasting negative energy densities can be sustained by quantum fields as 
shown by the example of the Casimir effect, modelling the plates as certain states of 
a full interacting theory. A computation along these lines was undertaken by Olum 
and Graham [68]; although there is a net positive energy density near the 'plates', 
their set-up maintains a negative energy density near the mid-point between them. 
However, it is possible that modified QEIs hold - see [T5J p. 176] for some discussion. 

• In terms of positive results, the averaged null energy condition is known to hold in 
general two-dimensional quantum field theories [H3]- In spacetime dimension of 3 or 
more, Bostelmann & Fewster [5] have proved that for a wide class of theories obeying 
the 'microscopic phase space condition' there are generally state-dependent QI type 
results on quantities that are 'classically positive', i.e., arise as the leading term in 
the OPE of a sum of squares. 

5.3 Singularity theorems 

I originally motivated the energy conditions by reference to the singularity theorems, in 
which they guarantee certain focussing behaviour. An important question is whether or 
not the QEI results provide sufficient control to guarantee that quantised matter also obeys 
singularity theorems. 
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Although this question is far from resolved, Fewster & Galloway [25] have recently 
shown that the hypotheses of the singularity theorems can be weakened to accommo- 
date bounds motivated by the QEIs. Unfortunately there is still a bit of a gap, because 
Hawking-style singularity theorems, concerning congruences of timelike geodesies, require 
the SEC (for which there is not a state-independent QEI) and Penrose-type results involve 
null geodesic congruences (which are not suitable for QEI bounds). Nonetheless, this is 
encouraging, and one can hope for more progress. For previous results along these lines 
see references in [2"5] . 
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